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ABSTRACT 


An analysis is made of in-plane and out-of-plane vibra- 
G2ons of a uniform filament which is bent at right angles 
and the ends of which are fixed. Shear deflection and rota- 
tory inertia are neglected. Curve sheets are presented from 
which one may find the first several frequencies of vibration 
in terms of the geometrical-.and physical variables. An 
amportant appiiecation is to the vibrations of piping systems. 

Also, as a separate project, a vibration program written 
by G. Fink for the CDC 1604 computer and modified by Y. S. 


Kim, has been rewritten for use with the IBM 360/67 computer. 
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£.’ INTRODUCTION 


A. GENERAL 

In order to-assure ‘long time integrity of nuclear piping, 
the applicable codes, ANSI B31.7 -— 1969 Code for Nuclear 
Power Piping and forthcoming edition of Section III of the 
ASME Boiler and Pressure Vessel Code, call for seismic 
analysis of such piping and require that the severity of the 
effects of seismically induced vibration be maintained below 
specified acceptable levels. Generally the problem of 
vibration analysis of the major piping systems is so complex 
that a major computational effort is required and this 
demands digital computer analysis using quite sophisticated 
programs. However, in the early stages of a design, when 
alternate layouts are being investigated, it is useful to be 
able to make some reasonable estimates of what a more elab- 
orate analysis would be likely to reveal. 

A recent document [3] concerned with assisting a designer 
to understand and comply with ANSI B31.7 contains material 
summarizing the vibration characteristics of several simple 
piping configurations. It suggests how the designer might 
employ this sort of information to assist him in estimating 
the vibration characteristics of more complicated configur- 
ations. This material should prove to be quite useful. 

However, the scope of the configurations about which 


such information is available is quite limited. One of the 


.s 


basic configurations with which a designer deals is a simple 
ell-shaped configuration. It was decided that it would be 
quite useful to determine and present in a concise form the 
Wignificant’ results relating: to the vibration of *such a 
configuration. Specifically, it is assumed that the two 
ends of the configuration are perfectly fixed or clamped, 
that the pipe is of uniform properties along its entire 
length and is of rather small diameter compared to its 
length (hence the word "filament" in the title of this 
thesis) so that shear deflection and rotatory inertia could 
be neglected, and that the right angle junction of the two 
legs of the ell is small compared to the overall dimensions 
wimeene ell Itself, so that 1G is not necessary to: consider 
the properties of the bend or elbow that would exist in 
practice at the junction. 

Such results are obtained and given herein. The basic 
parameter is the ratio of the length of the shorter leg to 
that of the longer leg; this parameter is designated by the 
Greek letter rho in what follows. A second important 
parameter, Greek beta in what follows, is required to relate 
the coupling that takes place between bending and torsion in 
the case of what are called out-of-plane vibrations. 

This simple ell-shaped configuration has the great 
virtue of being such that the vibration problem is naturally 
Givided into two disjoint sub-parts. One such part involves 
only displacements which lie in the plane of the ell; no 


worsion is involved in these cases. The other part involves 


both bending and torsion, but the motion of any point on the 
axis of the filament is only perpendicular bs the plane of 
the ell. This gives rise to the two terms used in what 
follows: "in-plane vibrations" and "out-of-plane 


wabrations." 


B. ASSUMPTIONS AND LIMITATIONS 

It is well to repeat here specifically the assumptions 

and limitations involved in this investigation. 

. The confticurdationm is ,formed from orteinalily straiehe 
uniform elastic members which meet at a right angle. 

2. The angle at which the two members meet remains a 
right angle throughout the deformation of the 
configuration. 

3. The ends of the configuration are assumed to be 
perfectly clamped. 

4, The effect of shear deformation has been neglected. 
Only bending and torsional compliances are con- 
sidered and deformations are assumed to be 'smaii. 

5. Rotatory inertia has been neglected. 

6. Poisson's ratio has been taken to be 0.3. 

Additionally, there are limitations implicit in the facv 

that values of the roots of the equations which are later 
Solved are obtained only to an accuracy of one part in one 
Million parts and in the fact that the pone ueee ene pro- 
grams used a mantissa having 56 bits. These limitations 


resulted in being able to obtain only about the lowest 


thirteen frequencies for a particular configuration. 


Also, 


only a few values of beta were investigated and there are 


some difficulties in the use of linear interpolation for 


other values of beta. 


©. NOTATION 


A, B, C,; D, F,H Constants in equation (2-5) and (2-14) 


E 


G 


&1> Eo» E3> Sy 


eee es ee ee ee et it 


a 


td 


or 


o- 


Modulus of elasticity 


Shear modulus of elasticity 


Functions defined in equation (2-6) 


Moment of inertia of cross section 


Polar moment of inertia 
Length of pipe 

Bending moment 

Defined in equation (2-4) 
Time 

Torque 

Shear force 

Distance from fixed end 
Lateral deflection 

Lateral deflection mode shape 
Slope 

Defined in equation (2-13) 
Defined in equation (2-29) 
Phase angle in equation (2-11) 
Angle of twist 


Angle of twist mode shape 


Constant in equation (2-2) 

Mass per unit length 

Mass moment of inertia per unit length 
Length ratio 


Circular frequency 


a8) 


Ii. THEORETICAL DEVELOPMENTS 


A. BENDING ANALYSIS 
Neglecting the effects of shear deflection and rotatory 
inertia the-equation governing lateral vibration of a uni- 


form beam is 


34 9° 
EI - a = 0 (221) 
ox ot 


mmere yt 1s mass pervunit length, EL is flexural rigidity, 
y is lateral deflection and t is time. Assuming harmonic 
eselilation with circular frequency w, the solution of 


equation (2-1) can be written in the following form: 
y(x.t) = Y(x) cosGat + wp) (222) 


m@nere j is a constant and.Y is.a function of x alone. sSub-— 


stitution of equation (2-2) into equation (2-1) yields 


4 
rar % ply (x) ai, oe oe 
x 
where 
2 
Yo pw 
P= ET (2-H) 


The general solution of equation (2-3) may be written as 
Y¥(x) = Ag, + Bg, + Ce, + Dey (2-5) 


where 


ab: 


SBICTAMA oi 


Sy 


(Cosh (Px) + Cos (=)! 


[Sinn ( Px) — Sin( Px] 


+ [(Cosh(Px) = Cos (Px) ] 


4% [Sinh(Px) + Sin(Px)] 


(a) 
(b) 
(2-6) 

Ce) 


(d) 


At any location x, deflection Yo slope Y » bending moment 


ML» and shear foree Ve Gan be expressed as 


EIP? (Ag, + Bg, 


A 3 


C(EIP?) ; 


+ 


Ag, + Bg. aa CE3 F DE} 
P(Ag. + BE. + Ce), + Dg, ) 


2 
E EP (Ag. + Bey, + Ce, + Dg5) 


Ce, + Dg) 


DP 


B(EIP?) 


(a) 
(b) 

(2-7) 
Ve) 
(d) 
(a,b) 

(2=8) 
(esd) 


Bxpressing the constants A; B, € and D in terms of Yo> 


1 


Yo> Mo and V 
! 
y; Pgs Ey 
M EIP“g.  EIPg 
3 2 
5) 2 
V EEP Ey EIP E3 


2 
E3/EIP 


&)/EIP 
Sy 


Pg. 
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9 and substituting into equation (2-7) yields: 


g,/EIP>| | ¥ 
B/EIP® 

(2-9) 
&),/7P M 


gE 
: 0 


The above equation expresses the deflection, slope, bending 
moment, and shear force at any location x in terms of those 


ac x = 0. 


B. TORSION ANALYSIS 
The equation of motion of a uniform cylinder undergoing 


torsional vibration is, 
—+~--—- >= 0 (2-10) 


where v is mass moment of inertia per unit length, J is 
polar moment of inertia of elastic resistance, and @ is 
Bue tee ot twista.. . Hor a. sinusoidal solution at cireular 


frequency w, 
Ox, E)MEe oC) COs CUt+*7) (2-11) 


where y is a constant and @ is a function of x alone. Sub- 


stitution in equation (2-10) gives 


2 
gtx) Ae Gay S10 (2a) 
ax 
where 
2 
ac = VO (2242) 
JG 


The general solution of equation (2-12) is 
¢(x) = F Cos(ax) + H Sin(ax) (2-14) 


and torque T is 


p= sgt) = JGa[-F Sin(ax) + H Cos(ax)] (2-15) 


ws) 


: Reto eon ah 
ry 4 ry 


zi ¢ Sas al 
aoa! . 
List) 


(63-45) 


teind) 


Coc: 


Para ae ; Ta H(JGa) (2-16a,b) 


Expressing constants F and H in terms of %o and To and 
substituting into equations (2-14) and (2-15), the angle of 
twist oe and torque a at any location x, can be expressed 


im terms.of T. and 9 as foilLows: 


0 
b Cos (ax) Sin(ax)/JGa}| | 


T -JGa Sin(ax) Cos (ax) ap 
0 
> 


C. IN-PLANE AND OUT-OF-PLANE VIBRATION ANALYSIS OF AN ELL- 
SHAPED UNIFORM FILAMENT WITH CLAMPED ENDS 


1. In-plane Vibrations 


In Figure 2-1, an ell-shaped configuration and the 


corresponding free body diagrams are shown. 


, Yi 
My 02 


Figure 2-1. Ell-shaped Configuration and Free 
Body Diagrams for In-plane Vibrations. 
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ae f= 
ae 


The ratio, (pi, Of jche dene 


th Of the shorter Jee to. that 


of the longer is a basic geometrical parameter in the 


following analysis. 
At the clamped ends defle 
1 t 
mi} 02 
portions, deflections Yh, 


4 are equal to zero. 


1 and 


cation of equation (2-9) to t 


B ES 
grips, 0+. prp3 9+ 
Ey & 3 
— /M + V = 
erp 91 ype 
Ey 
toi eke) On) 
g g 
3, Le = Vga > 
EIP EIP 
: Zz 
~ Moo t 5 Nga 
EIP EIP 
6. Mase as Veo = i 
we eee Le 


etions You? Yoo and slopes 


Also at the junction of the two 


Yh are equal to zero. Appli- 
he two portions yields, 
0 (a) 
ie (b) 
L1 
Ce) 
(2-18) 
0 (a) 
te (e) 
L2 . 
® 


where a indicates evaluation of g, (x) for argument (PL) and 


Bs indicates evaluation for argument (pPL). 


At the junction, 


¥ = Y 


' 
L2 


(2-19a,b) 


Combining equations (2-19) with equations (2-18), the 


following set of equations is obtained. 


15 


g3/EIP® 0 g5/EIP? 0 


0 z4/EIP* 0 Z,/EIP> 
% a) ; =0 (2-20) 
g,/EIP -&,/EIP &2/EIP -&3/EIP 


For a nontrivial solution, the determinant, D, of the above 
matrix must vanish. Elementary manipulations on D lead to 


the following form. 


20 (2221) 


The terms of equation (2-21) are functions of op. For 
any value of p, the zeros (PL), of the determinant can be 
found. The corresponding frequencies are evaluated by using 
equation (2-4) as 


EL 


2 
Gi Oe ioe 
i ni Tew 


(2-22) 


Solution of equation (2-21) will be explained in 
Section III-A. 


2. Out-of-plane Vibrations 


Out-of-plane vibrations can be analyzed similarly. 
However in this case the effect of torsion must be taken 


into account. 
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In Figure 2-2, free body diagrams of an ell-shaped con- 


figuration for out-of-plane vibrations are shown. 
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Figure 2-2. Free Body Diagrams of an Ell-shaped 
Configuration for Out-of-plane Vibrations. 


At the clamped ends deflections Y and slopes 


G22 92 
t ! ° 
Yo. Yoo and angles of twist %o1> do0 are equal to zero. 


Poplication of equation (2-9) to the two portions yields 


aly 


a) 


t¢ 2 


g g 
223 yy + 2 V (a) 


Ald: ose alee aas 

Yr = ot ton — Voa (b) 

DE Vs oe Gacas “i von (ce) 

fines Soe son (d) 

me 2a + Eo Vv (e) = 
Op REE AU aa PEEPS OE 

ae nef yp 3, Vigp = 

Mio = 81Mgo + a Yoo (g) 

Vio = P&Mos + 81Vg9 (h) 


where g, indicates evaluation of g, (x) for argument (PL) and 
g; indicates evaluation for argument (pPL). 


For the torsion, equation (2-17) applies and gives, 


. Sintal) 
(eg a ee SO (a) 
aa ale 01 
T1* Cos (aL) To1 (b) 
(220) 
Y-sin(aps) 
on = Ge ann OP (c) 
i> = Cos(apL)T), | Cd) 


At the’ junction: 


18 


SGaB601T) Irmeey: 


M =o : 7 eat €acb ) 


Eee or Ll L2 

Sos : _ vt _ 
Pra) pee sn he (c,d) (2-25) 
Md, eke Woe ae eee (e,f) 


specializing to the case of circular cross section and 


Sosuming Poisson's ratio is”0.3, ‘then 


a es eS (2-26) 


Combination of equations (2-23), (2-24), (2-25) and 


(2-26) givesthe following set of equations. 


= 1QP oe 
0 yP ~<~—Sin(L) 0 0 
Les 
g,P 0 0 ~s—Sin(apL) E3 
ss = 0 
Ee &5P 0 ° 8) 
Ey 
Ey 0 O- -Cos(apL) = 
0 gE, Cos (aL) 0 0 
C2227) 


For a nontrivial solution, the determinant D of this matrix 
must vanish. From equations (2-4), (2-13) and (2-26) it is 


found that 


(aly = 6(PL)* (2-28) 


V9 


, 
¢ 


ei 
. 


a 


ay-¢ ) 
rin gi-n 


c 
‘ 


> 
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a] 


where 


B = = (2-29) 


is a dimensionless parameter. Substituting equation (2-28) 
into equation (2-27) and performing elementary manipulations 


on D give the following equation 


3-63 E : Bo 8; 
ee, -= 3 sin[6(PL)*] 0 ore. 
E, 0 ) 223 sin [p8 (PL) “1 She OT 8 0 
fe ..& 0 0 GB, By 
(2-30) 
g, 0 0 (PL)Cos[Bp(PL)] g, 0 
Meee, (PL)Cosis(PL) } 0 0 By 


The terms of equation (2-30) ren tinetions or 6 and 6. 
For given values of 0 and 8, the zeros (PL), of determinant 
D can be found. The corresponding Prequeneies are evaluated 
by using equation (2-22). 

Solution of equation (2-30) will be explained in 


Section III-A. 
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TIT. PRESENTATION OF RESULTS 


A. METHODS OF SOLUTION 

The equations (2-21) and (2-30), developed in Chapter II, 
can be solved either by algebraic manipulations or by using 
direct computer solution. 

ina first approach to a solution, algebraic manipula— 
tions were performed to reduce the order of determinants. 
Then transcendental equations were obtained and two computer 
programs, one for the in-plane solution and one for the out- 
of-plane solution, were developed to get the roots of these 
equations. 

In the computer programs the transcendental equations 
were used as function subprograms and roots located in an 
interval. Then the Newton-Raphson method was used to obtain 
good approximate values of the roots. 

The procedure explained above led to some difficulties. 
During the solution of the equation which was derived from 
equation (2-30), the Newton-Raphson method failed to give an 
expected root because the slope of the transcendental function 
reached very high values and this caused a jump to adjacent 
roots. In one case, the number of iterations required 


8 (for 


6 


exceeded three hundred with an allowable error of 107 
the quantity PL) and increasing the allowable error to 10. 
did not make much difference. Finally, comparison of 


results with the classical solution in the case op is 


Bok 


“4 ie 
‘ 
7 7 val 
— — > 
7 edhe - 
-‘S - a 7 
ans 
nt 
— 7 Be as 
ve 
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oe 12h eaarto ef} utd tw 


equal to zero showed that some of the roots were 
extraneous. 

The difficulties mentioned in the preceding paragraph 
‘were eliminated by making direct computer evaluation of 
the determinants (2-21) and (2-30). For this purpose 
subroutine BRODE was developed and each element of the 
determinants was stored in this subprogram.. In addition to 
subroutine BRODE, library subroutine DTERM which was 
developed by Jean Bow in 1967 was used to evaluate the 
determinants. This modified procedure resulted in reduction 
of computer time required and in the elimination of extra- 
neous roots. 

Difficulties with occasional failure of the Newton- 
Raphson method were reduced by using the method of successive 
bisectioning. In this method the number of iterations 
required for one root did not exceed fifty with an allowable 
error of Tome. This procedure also resulted in conservation 
of computer time. 

Detailed instructions for use of the computer programs 


are given in Appendix B. 


B.. (GRAPHS 

Numerical results given to a large number of significant 
figures were obtained primarily for the purpose of investi- 
Bating the accuracy and integrity of the solutions. Chapter 
IV will give a discussion of these investigations. 

However, so as to present the results in the most con- 


wenvent. form for practical use, graphical outputs 
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were also obtained by use of the library subroutine 
DRAW. 

One graph was obtained for the in-plane solution. For 
the out-of-plane solution five graphs were obtained cor- 
responding to five different values of 8. In Seah ease, 
the ordinates are values of (PL) and abscissas are values 
of p. On these graphs the first (lowest) curve is for the 
first mode, the second curve for the second mode and so 


on. Graphs are given in Figures 3-1 through 3-6. 


C. SAMPLE PROBLEM 
A sample problem will be solved to show how to use the 
erapns given in the preceding article. 


Properties of ell-shaped configuration: 


Ly aa eS 

L, 2 PLO l= oe" 

Outside diameter of pipe : 6.625" 

Inside diameter of pipe 5.065" 

Wall thickness of pipe Ose s0" 

Weight of pipe s 18690 Ib/re 

Thickness of insulation 20 Rae 276" 

Weight of insulation ote. 902 1b /if se 

Weight of contents 112351) 1B/ee 

Medulus of elasticity E- :°° 27,400,000 1b/in® 
G : 10,600,000 1b/in@ 


For the use of graphs, values of p and 8 must be known. 


Por this problem, 


25 


(adh) JO SOnTeA 
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In-plane Vibrations. 


Figure 3-1. 
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Figure 3-2. Out-of-plane Vibrations, Beta = 0.00001. 
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Out-of-plane Vibrations, 0..005.. 


Figure 3-4. 
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Figure 3-5. Out—ol—-plane Vibrations, 6. =-¢.0n. 
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Values of Rho 
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Out-of-plane Vibrations, 


Figure 3-6. 


ae, 


(3-1) 


where v is the total mass moment of inertia per unit length 
fer the torsion and } is the total mass per unit length 
undergoing lateral vibrations. Assuming that the contents 


(water) do not rotate, then 


= ; + WV, ‘ 
. Ypipe Vinsulation 


Ueotal — Uoipe i. Hinsulation * Ucontents 
numerical values of Venton and He Gea are 
as 2 
Vegetal 0.0509089405 lb-sec 
a 2 ee 
He eal. 7 1.068875490 lb-sec*°-in 


and from equation (3-1) 


B = 0.0150806896 


For evaluation of frequencies, equation (2-22) should be 


used as 
w, = ay EI 
uL 
i 
then 
Ba. ye -1 
eee Tis 8.221317961 sec 
WL, 


30, 


as 


{ fi 
1s y 
} Se 
i) VL 


—j 


its 5 


wy 


To solve the problem, values of (PL), can be found from 
the graphs. But no graph is available for the value of 8 
found above. Therefore interpolation is required between 
the values taken from Figure 3-5 and from Figure 3-6. Values 
of (PL), are given in the table below where the exact 
solution was obtained by a computer solution using the 


specified value for 8. 


MODE EXACT SOLUTION OBTAINED BY 
NO. SOLUTION LINEAR INTERPOLATION 
i) 2.624579 2.624302 
2 4.533513 4.521991 
3 6.729044 6.706905 
4 7.862686 817582 
5 10.518163 10. 230722 
6 11.9067 33 11.698410 
7 13-1 50277 13.442683 
8 14.259158 15.619057 
9 16. 376500 16.379861 
10 17.185014 17.403866 


As seen from the table, for the lower modes, linear inter- 
polation gives good results. However there is greater error in using 
linear interpolation for the higher frequencies. One can 
make what appears to be the reasonable conjecture that if 
one were to plot curves with (PL) as ordinate and 8 as 


abscissa and with different curve sheets corresponding to 


Si 


. oe 
Ls 
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fixed values of p, there might be some crossing of individual 
eurves; that is ,. for) certain: values of 8B ere Gig boat modes 
might have the same frequency. Clearly this would confuse 
the significance of characterizing different modes simply 
by ordinal number in the order of increasing frequency and 
would indicate that a more sophisticated interpolation is 
required. However, no such graphs were prepared and there 
do not seem to be enough values of 8 represented in the 
Braphsvat hand’ to permit making repilots.” Thus’it is 
suggested that some future investigator might find it 
interesting to employ the computer programs, incorporating 
slight modifications, to obtain such graphs and thus inves-— 
tigate the reasons for inaccurate interpolation for the 


higher modes. 
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IV. ACCURACY AND ASSURANCE OF SOLUTION 


The allowable error permitted in determining the roots 


6 


and this was assumed to be satisfactory for practical 


The allowable error can be changed from 107° to 


510° 
purposes. 
any desired value by changing the value EPS in the computer 
programs. This can be done by changing the card 000090 intte 
program IN-PLANE and card 000100 in the propran OUT-OF-—PLANE. 

Assurance of solution was established in the following ways. 
1. In the case p is equal to zero, the ell-shaped con- 
figuration becomes a simple built-in beam and solution for 


this case can be compared with the classical solution. 


4.7300407 
7.8532046 
10.9956078 


7300289 
8531849 
-9955803 


73004 EL 
7.8532047 
10.9956081 


14 


A 
2 
5 
4 
5 
6 
7 
8 
9 


= 
>) 


co 
Obtained by computer solution 


» 1371655 
aye 
210% 
23. 
Ge 
29. 
32. 


2787597 
4203522 
5619449 
7035376 
8451302 
9867229 


SOs (PL) Cos (PL) = 1. 


pia 71302 
2787165 
4203012 
.5618858 
7034814 
. 8453738 
.9890104 


14 
L¢ 
20 
23 
26 
29 
Se 


of the equation 


«1371655 
~2787594 
4203526 
5619449 
1035412 
-8451861 
9866973 
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In equation (2-21) when p is set to zero the value of 
the determinant is zero regardless of the value of (PL). 
Therefore op was taken to be 0.00001. 

In out-of-plane case for the larger values of 8, pure 
torsional modes appeared and this can be seen from Figures 
(3-5) and (3-6). The values of (PL) given in the above 
table were taken when 8 is equal to 0.001. 


2. Consider the casesshown in Figure 4-1. 


(a) (b ) 


Figure 4-1. Cases Having Reciprocal Values of p. 


The configurations shown in the above figure are identical, 
therefore two systems must have the same frequencies. Using 


equation (2-22) as 


(PLYeI (Er (Pee = 
ae a 


then 
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«3155336 
2490132 


- 7635958 
- 707676 
~540341 
-024013 
-987974 
. 8231918 
- 3082077 
2712829 


- 9832348 
- 7014996 
- 3772076 
- 8416393 
- 8159008 
- 6035633 
-0968311 
0625606 


. 8673419 


- 35007 03 


. 1877668 
- 6245066 
. 3817979 
23530301 
~7701705 
.5120064 
- 4939872 
91 15959 
-6541038 
.6356415 


For out-of-plane solution in the case 8 


1.4916187 
2.300753 
3.6886158 
- 4208958 
-4080044 
8018597 
-5488148 
.5315098 
~9339245 
-6796268 
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-000000000 
-000000000 
- 000000000 
-999999999 
-000000001 
999999999 
-000000000 
-000000000 
-000000001 


-000000001 


L- 999998209 
-999996916 
-99999 3485 
-999965571 
-999980037 
-999977062 

-999894233 
.999946210 
-999948965 
9991 TT 515 


= 
= — = 


3. In the case p is equal to one, it is expected that 
among the in-plane solutions will be many for which the 
moment at the joint vanishes. Thus, these cases also per- 
tain to the clamped-pinned simple beam. Now two clamped- 
pinned simple beams of different lengths, & and 2£*, may have 
matching frequencies in their spectra. That is, for example, 
the second frequency for one may equal the first frequency 
for the other. Specifically the following case was 
considered. 

The third frequency for an equal-leg ell is equal to 


the second frequency of a clamped-pinned beam of length &, 


namely 
2 EL 
OCP 
s) 3 wer 
where 


(PL) 3 = 7.06858275 


The first frequency of an equal-leg ell of length g2* is 
equal to the first frequency of clamped-pinned beam of 


length 2*, namely 


= 2 EE 
Piers eee 


where 


(PL), = 3.92660231 


mhe ratio 0 was determined so that then the frequencies are 


equal, thus 
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Then performing an analysis for an ell-shaped configuration 


with the value of op, found above, the second frequency is 


indeed, 
EL 
Wo = (PL), [—F 
ys 
where 


(PL), = 7.068384686 
This value shows close agreement with the value 
(PL). = 7.0685275 


given above. 
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APPENDIX A 
DESCRIPTION OF PROGRAMS 


A. GENERAL REMARKS 

Programs IN-PLANE and OUT-OF-PLANE are FORTRAN-IV 
language digital computer programs designed for use in the 
vibration analysis of an ell-shaped uniform filament with 
clamped ends. 

Both programs employ double precision arithmetic. Sub- 
Foutines are provided in binary form. Trigometric and 
hyperbolic functions, available in the IBM-360/67 library, 


were used. 


B. PROGRAM STRUCTURE 

Bach program consists of a four section main program and 
tao subroutines. The function of each section of the main 
program and of each subroutine is given below. 


1. Main program subsections and their functions. 


INPUT Controls read-in of data and allocates 
storage locations for arrays and subscripted 
variables. The allowable error for deter- 
Minavion Of roous. 1s specified in this 
section. 


ITERATION im thise section the roots. of the determinant 


AND are located in an interval by searching the 
EVALUATION sign changes of the determinant and then the 
OF ROOTS method of successive bisectioning is used to 


evaluate the roots. The maximum number of 
iterations was limited to three hundred and 
this was assumed sufficient. If the number 
of iterations becomes as large as three 
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OUTPUT 


GRAPH 


hundred the process terminates and the 
approximate root is tagged to indicate a 
possible inaccurate result. 


Controls output format. 
Controls the graph output and uses the DRAW 


subroutine provided in the IBM-—360/67 
Library. 


eae oubroutines and their «functions. 


BRODE 


DTERM 


In program IN-PLANE equation (2-21) is stored 
and in program OUT-OF-PLANE equation (2-30) 
is stored. Evaluates each element of the 
determinant and obtains the value of the 
determinant by calling subroutine DTERM. 


This subroutine computes the determinant of 
a matrix of real numbers by Gauss! method 
of elimination with row pivoting. 


PROGRAM NOMENCLATURE 


R3, R4 


Elements of determinants. 

Checks sign of determinant and if there is a 
sign change moves control to the method of 
successive bisectioning. 

Values of 6 ian subroutine BRODE. 

Checks the acceptability of roots. 

Value of 6. 

Single precision version of RO. 

Value of determinant . 

Allowable error in roots. 

Value of determinant’. 

Index of modes (Ordinal number). 


Values of (PL) used to search interval. 


Values of determinant in the method. of 
successive bisectioning. 


Values of 0. 
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ROOT 


XX2 


Roots of determinant (PL). 
Corresponds to RO in subroutine BRODE. 
Single precision version of ROOT. 


Middle of interval. 


Upper and lower limits of interval contain- 
ing Che root. 


Values of P, xX, XX1 and KxX2 in subroutine 
BRODE. 
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APPENDIX B 


INSTRUCTIONS FOR USE OF PROGRAMS 


A. DATA REQUIRED 


ie 


Data required for program IN=PLANE 
a. Number of p's (NRO). 


b. Number of modes sought (MODE). 
ec. Graph control (NGRAPH). 

d. *<Values (offi. 91 (RO). 

ebe Graph title, CETITLE): 

Data required for program OUT-OF-PLANE 
a. Number of o's (NRO) 

b. Number of @"s (NBETA). 

c. Number of modes sought (MODE). 
ad. Graph control (NGRAPH). 

e. Values of 6 CBETA). 

f. Values of p CRO). 


g. Graph title (ITITLE). 


B. DATA DECK ASSEMBLY 


Li 


Program IN-PLANE 

The data deck is ordered as follows: 
a. Problem statement card. 

b. Valves of 0 card(s). 

¢. Graph tavie (cards. 


Program OUT-OF-PLANE 


The data deck is3erdered as follows: 
a. Problem statement card. 

Bow Valueseot 6 card(s): 

¢. Waluesnofiep: card(s); 


ad. Graphteitie cards: 


C. LIMITING THE NUMBER OF ITERATIONS 

In the programs the maximum number of iterations was set 
two 300. If the number of iterations reaches 300 for any 
root, the current value is assumed as a solution but in the 
Sucput it is marked» as ITER. EXCEEDED 300. A user either 
accepts the root or changes the cards 000430 in program 
IN-PLANE and 000460 in program OUT-OF-PLANE by changing 300 


to any desired value. 


De. GRAPH 

A user can get graphical output by option. DRAW sub- 
routine, available in IBM-360/67 library was used and three 
CALL DRAW statements control first, intermediate and last 
curves to be plotted on one graph. DRAW subroutine requires 
Single precision arithmetic. Therefore variables CEM and 
SEV are introduced to change the values of RO and ROOT from 


Souble precision to single precision. 


E. FORMAT OF DATA 


1. Program IN-PLANE 
a. Problem statement card. 
The Gard has 3 fields under the control of 


following field specification. 


Ho 


q 
% 


Sata 

Field names and the information each field conveys are 
given below in order. 

Field names and content. 

(1) NRO Number of p's 

(2) MODE Number of modes sought. 

(3) NGRAPH Code specifying whether there is 

tor be graph output. 


(NGRAPH O "No graph output) 
(NGRAPH 1 Graph will be output) 


ood 


bu - Values of 9 card(s). 

Field speciticatrion  7FiL0'.0'. 

If there are more than ( @gittrerent values of 
op, additional card(s) will be required with the same field 
specification. 

eC. Graph ticle cards. 

The title of graph is read as data with the 
field specification 648. For the title two cards are 
Zmeogurred. The first card corresponds to the first line of 
the title and the second card corresponds to the second 
fine. 

2. Program OUT-OF-PLANE 
a. Problem statement card. 

The card has 4 fields under the following field 

specification 
ho 3 


Field names are given.below in order. 
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(2). NRG Number of p's. 

(2) NBETA Number of 8's. 

(3) MODE Number of modes sought. 

(4) NGRAPH Same as in program IN-PLANE 
Values..of 8 cards. 

Field specification 7F10.0 and field name is BETA: 
Values of p card(s): 

Same as in program IN-PLANE. 

Graph title cards. 


Same as in program IN-PLANE. 


Ky 
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APPENDIX C 
FLOW DIAGRAMS 


A. GENERAL 
Flow diagram of the program OUT-OF-PLANE is given. This 
flow diagram includes all parts of the main program. Flow 
diagrams of subroutines BRODE and DTERM are not given and 
their functions were expldined in Section A-2 in Appendix A. 
Program IN-PLANE is similar to OUT-OF-PLANE, therefore 


its flow diagram is not given. 


B. FLOW DIAGRAM OF PROGRAM OUT-OF-PLANE 


Flow diagram follows. 
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READ: NRO,NBETA 
MODE, NGRAPH 


READ: BETA(IT) 
ih. = i, NBETA 


READE ROM SS). 
Jue tale NRO 


CALL BRODE(P(1), 
RO(J), 
BETA CGE), DET) 


F(1) = DET 


CALL BRODE(P(K) 
RO(J), 
BETACE > DET) 
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R¢K) = DEP 


B =F(K)*F(K-1) 


0) 
oO 
(S) 


Le 


KG) St C1) 2 
+(XK?2 72 


BBB = XX2 - XXX1 


< 
=f —-ROOT OR) 


‘i 


CALL BRODE 
(Ax TROCS ) 


BETA(1) ,DET) 


Rs = DE 


CALL BRODE 
(X(T) ,RO(J), 
BETA(L)DET) 


pe 


on 


Se 


WRITE: ROOT 
(NJ By. TPER: 


EXCEEDED 300 


eat 
5 | 
<Eenmet 
P(Ktl) = 
PECK) Soros 
<n 
ge koe eh 
WRITE: , 
ROOWON i 5) 
1 ,MODE 
1 ,NRD 
1,NBETA 
NGRAPH:0 8 


f 
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READ : ITITLE(T) 
i aka i 
READ: ITITLE(T) 
ae ise 


MM = 1 
MMM = l 


sa - 
ROOT (M,MMM , MM) 


<i n> MMM = MMM + 1 
<a mci CALL DRAW 


CALL DRAW 


CALL DRAW 
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APPENDIX D 


PROGRAM LISTING 


A. INDEX 


ice 


2. 


Program IN-PLANE . 
Main program: 
Section 
INPUT 
ITERATION AND 
EVALUATION 
OF ROOTS 
OUTPUT 


GRAPH 


Subroutines: 
BRODE 


DTERM 


Program OUT-OF-PLANE 
Main program: 
Section 
INPUT 
ITERATION AND 
EVALUATION 
OF ROOTS 
OUTPUT 


GRAPH 


Sequence number 


000090 


000170 


000590 
000810 


001060 
001310 


Sequence number 


000100 


000190 


000630 


000870 


pie 


- 000130 


= 000550 


- 000770 
- 001030 


- 001300 


- 001610 


- 000150 


- 000590 


- 000830 


- 001110 


2)) 
55 


58 
58 


WP 


Subroutines: 


BRODE 


DTERM 


B. LISTING 


001140 - 001600 


001610 - 001910 
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PROGRAM IN-PLANE 
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PROGRAM QUT-OF—PLANE 


eWlolelelsjlelelelololelolojlolalolejelelolololcelejlolalelelelelelolelolelelolejelelelelelel2) 
FANMOFINDOE DRHOANMTNOEDMNOBNMTNOEDHOGANMSFNOEDAORNM TINO 
DO 000 CO OO RSS AAS RRNA NAN NIN NN AI 909.09) OH OOO OO SS St 
lelejoleleloelelolojolelojolejolejololclojololejelolelolojoelololejelojelelelelelejejejelejeye) 
eloleleleleloloeleloleloleelolel el el elelelel cel eoelelololololelolelelelolelolelelelelelele) 
lolol elelelelelololelelelololelelejele} ejlajolelojeleloalclelelelelelelolejlolelelelolelelelele) 


)+X(350),RO(60) ,ROOT(15,60,5),BETA(5) »SEV(60) »,CEM( 


a0 Ao a A NO) 
4 


INPUT 


OOO 


ee 
woo 
Or oe 
le te @) 
Oo ™ at 
, ota el 
Ost NuN™ 
oO ow & ow 
ek LIN) 
| pa Yael <a ea Ye 
NaIedade 
WoL aa 


RO»NBETA,MODE »NGRAPH 


et N 


oc 

mt 

(Sys) 

edad 

— uw wi he WwW 

‘) QO (=) Ww Oo Ss 

oO . bas (=) CI SKS) 

w ~ ~ oe) om PETE) 

—J =) =~ + OO 

Ww — —_— = | — eo°o 

(se) <— — ~~ <td oOo 

ke ke ater ee 

Zz uw uw Re W Ke 

CO a m (o) NW © OO 

— « WN fee) - ee 
= = =~ eye — tt + 
<a =; az) Oo (Sy (BS ele ne s(os ex 

= = ~ te Qe =e 0 
- OQ O oO = © ae (=) 
<I (og aw (S} e-O0 % 22 ome 

> - . Oo NAOw ey << 
Lu og ~ - NX _-— °° Oo 
re wi 4 ~ ~<a et OO o 

areow = o~  — CS NIOxX ~ AO 
Fe (OE(aa) Qa Noa AHO <nNax «KX OO -~ 
<< 2a Om O~— io) <>< ee (S} 
ee OW OCww Yo ~~ +KM ew WI OO (‘S)} 
Z eae WO ANNO ~O A wFimwlUO A oe oa) 
Cc nin COrFrRNOrFWL o 1—=x<xNIO O FRO © 
(al eee!) eC et eu e YYKM eY CLC Odes GD 
OMCO DAO ~we~wKIOOrDOF © Owe wl 
<< Oxco) ih ie SiS eo -caee Welco StH) Swimm ee< “se 
CO wae Kderendtse~O Nm HOYOHOC AIRE Ne 
uu Aer N Slee tr Ow HON OR 
- OCOONn~s~~OM~ fl Lh WK KS OU IMIS ULUKO-IL 
—= OOZAOLAOOW Osx kK KOMOCOUC RS KOx 
Oo On 
a 


OVUO 


Dt 


~ 


elolelololelololelolelelole) elo) elololelololeloleoleloleloleleleleleleleleleole) el oleleole\e\ela\e) 
MSODOANMAPOIOHDHAOANMASFHNOE DHOANMFNOHDHOAUMFIANOE DHOOM + 
PFTFIOOMNMIMINWDMADIN DODODDDDODDOEEE EE EEE ODDDDDDDDDOAHAHAADA 
eleleolololelolelolelolelolejelclololololol Soelelololololelslelelelelelelelolelelelelelelelole) 
eleleleleleleleloleloleleleloelolejlelelelolelelel Sele] elelel ele) eee) o] ee) 2) eel o) eel eee) 
elolelelelelel el elolol eel Plelelelolol eeolelelejlal ele) el ele] elolel Sol eee! 2) elelelelelelelelo) 


~ 
oO 
(2) 
m 
oO 
WwW 
(=) 
uw 
WwW 
O 
x< 
uj 
e 
[ad 
WwW) 
= 
— 
ae 
es) 
te | 
2) 
x 
e 
—_—O 
ao 
eu. 
=2 =O 
ex<re 
Zo 
--0O 
Ee © 
ON (=) 
mM Hee O 
~ ~— XxX e 
x xZR=~ Oo 
W i ww + 
_—~ -—— eQD — 
J Jno vv 
bad AO=T ~~ 
"O07 = e ALLW 
Dent “0~G IlDDD 
2 awe axZzZ2Zz 
BOW Om WT erie 
FERRERS Z+ +E RK 
— O ONe~Zzvzzz 
VwOOOOCOUWwI~ODO 
SOOO FLSsZAOUU 
N YT WO OFOD 
et aoe Nee 


OUTPUT 


OOO 


) 


+ 
* 
% 
* 
H* 
HH 
+ 
#* 
HH 
HH 
% Ht 
+ * 
% #% 
— itu 
oO #*# ~ 
tise * 
QD Kt * 
a #+ + 
H 3 * 
% 3 + 
ce ## * 
Ww + * 
(oe) Gas + 
= *++# + 
wHtte = — + 
O Fx WM Ww + 
wu £t#Ht 2 e + 
Q **x O - % 
Se cap bey “Ww + 
t+ Het & ~ 
4 +e « i 4% 
A cps [ee —a— + 
<x #H © ee * 
Oo +e Uy + 
aA tH > oO + 
- #t a5 pS 
a a: S| iy te 
Cal sy es ~ *€#+ Aa -—— 
ce #H Ft <x + TC wo 
Ww #Ht A WNW + N eet 
oO tH tt a — en tt tad ~e 
HH SN < 0% 
HH LW —e * O wR 
=z *# O —~ + es ~LL 
— 6 aIut + e Ze 
tH FH k- hm ~ HH WwW MX 
LOS wil +t Aa ~o 
oO tt O ao *+ QO Fea 
a tte ds fieeY., a= Cie 
N 3 ee) EWN 
a Ht NOWIM Mm TOCR 
~ T+ OY Oe ee “OO ee 
< Ox aw << OK SS Oe 
Wy ett ae eC OO OO aomf- 
ST AH WOrtneresty ts Ta + 
Sal ot ~~ il « o eine 


awe ew tee Dime mew me ZMe 
MAtOKMOORHHOOH- DORK HO 

woo ere tt ees co Se fe ee Will 
0O~ O~ HF OM COU~0O~0~O0O~DDD 
wee HF OO em Ew rr OK Zee 
Wott HH WINARW KD I TOU 
RE rE Hr SOW RSE ENR Sere 
Ve 4 See ® Sr deel Cho Aen © Ae] ° SU | 8 
COCO COOOLOCOCOOCOO0OO 
SUL SLL H# SILOOSILSWSLOSLOOU 

et 

mm + uw OO = & MOOO 
AOOO 

Onn 

LOY LVL 


58 


GRAPH 


OOO 


MM.MM) 


ees 


WOO- ll =O 
eOOO | =STe 
Do wodtl=z= il 
aw two - 
Ine ewOO00O= 
LE OCO= 
O~mwn SIMONE 
ZO0O=MMOM~— 
~adae = 
WwwWOoOOOW 
Sse eLoOooOoWM 


° 
oO 
. 2) 
wT 


io@ 


w 


M 


wa 


lelolololololelelole] eee) ejelelelole@) 
WOM DHORANMTFODOEDMDORANMY 
AAARDOOODODCOC C00 ASS 
feolelololer te ee _ 
DOOCVCCCOOOOCCCOOCO0O0O00 
elelelolololelololelejlele]ejlele]elele) 


rs rs 
ol! o 4 
Lod tal A ad 
fo) [e) o 
et c=! ed 
tad A od e 
= ~- ~ 
oO oO o 
cS Cs Cy 
oO oO oO 
tad td tol 
oO oO oO 
o 
oO (S) oO 
e o 
oO (2) oO 
e e e 
wh ita) iva) 
Ld e e 
N N N 
e e e 
oo oO oO 
Oe - 
mu Ww uJ 
gee | —) =) 
- ke 
On mm a 
ee ke Ld 
oO tem i — 
Oo Oe « « 
at O — —) 
+ Ww uy uw 
—co (ea) faa) 
oO aa <a o = 
| a oe =!) —) 
We Co ~ 
(SS) leek) oO (o) -~ 
Oow e ~ N 
oO es N ina) _ 
NO & Co rs 
atu> > > = 
e el uw Ww rT] 
COzNH ~Y Y - 
Wes e—  & - ew 
° = = ae I< 
= rEtoyeyte Ww Ww dao 
= ZzOZO O O Hse 
= e dae CY ~ ° 
= A= OQ oO Oe x< 
~ awwe fag oO VOU 
SS <aooz Zz =z *Oea 
as eO0w~ Ow Ow WO 
mse =O HOS es 
waWeeecat tot Wao ert LU LU 
SDOwor te Fore O~DD 
SZWZWUO an) QO wke2Za 


Samecee K=O SO alice 
ed ae | ad ee | ad ee | od ed ee 

SZ~ee~ssN IN UHR ZZO0O 
WOWUW WI AoOatOdadewOOOrZ 
OOK RRSOIUOOIJOUIUSLOOMNW 


ei ee ee 
o oO ° Oo OOCO0C°O 
Oo 2) (=) oOo O9O000O 
- a N MY TFUINMO 
na) t+ + Fr TSO 


eloleleloleleojeleloleleloelelolejolelolelelelelelele) 
FNOM™ DHOANUOPFNOEDRHROANMTNOF ODO 
Se tt RSI NUINI NNIAI NI N ANIA MM MMH Meo OH + 
ren Koen Knee Aen aren Kaen Aon Kae Lame hones fone kamen Lane Looe Ramee Laren tee Lee Raves Lae Raee Kaen han here Ree | 
GDODCOOCOOCOCOCOOCOCOCOOCCOCOCOOOCOO0CoO 
Wolo eelel 2lolelelelelelel el al elelelelelo/ el eelole) 


ed = Saad -_ 

Ww) ~_~-_ = = 

~ Ss >> = 5. 

(oa) ad —~ DD ~— 

am NY Pare 00 = 

en (oie) Sa OO _ ea) 

a | OU OW meN 7a) faa) 

ae ee QAO NO ro) ~ 

oie 1s tt (eye) + oO 

>t -—— a~—a tM _ = 

>I ~D ~D. — e > to3) 

~<a > adil -S= De > e 

QW oO mM BZ a * 

OX & NOD Ze 2s az ~ 

ode (Sia) NA Ma oY) oe en 
wa QOw~ Qw~ Qr “AQ ae oo 


N 
St ae * Ko %+>- +O a+ ow ~ 
e FINNS INS INN I ENNOSNGS 
Re OOY > Hot eT a a a a aaa 
OMZNOHKH KH CANMTONOANMSOOAINM TNO 
Ot © o> Ot we hr RR wR RR OR RO RHE HR HHH 


OMHOyYFNDS>addadadddaqdaqadqadqaqdaqddqdaacd 


a9 


se 


Yin s — aor rt ps) ATS p 
e -— _ 3 
ef eOE eT yes Dy e2eS.0,3 7 ae BAD. eT te 
: (TZA)i 
: 7 OHHe ro 
# EOE aT sO OG oO uO ote Se0y FITIT2  SIBAS 06S. VIZ .MaD OND were 3 
ef OL eT per BeOS Petr SeOe STITT y J98451.0, 6, WRe,4) 
(S¥e*<TdA 1). ke ccHl bs 
W008, 88,YVISTGRS 3 
;* =4 7 | 
6 - ‘ 
Levy) 4090-1 vv 14 t PP 
$ipcov0~{ Ue Oi ber- 
* 47 Paaqr 
Ti PR ieiHt vy pried seve pe eh) 
Hin Fe LOK CtlpeAIeOler—<ts, 1} 14 
Let vs 7A 
TU Cee ee ee / 
: BUSoe sat liattonc it?s! 
a > id. ape 
i i. et he Et 
8 Dy 
DL te ait 
i: a¢@ 
_ 7 At atte? ‘ 
* 
7 erty Hy: bia . 
: ‘ 


_ 


COCCCOCOCOCCOOCOOCOCOCOCO0O0o 
ANMTNOEY DHDOANMSTNOYDDHDO 
PSSST SFSTMOMMNOIMNNAMNININO 
Ssqqadt aaa ss eS et NA eae et Se 
eleleloelelelelolejlolelele)2ele)elejole) 
leoleloloel elelelelolelelole)lejeleleleleye) 


H 
H 
COSCV) 
A,DET,6) 


Or NO 
OQ e« % an) 
—~vTt ~m tT 

== ~—o 


~— —> 
1INN>><tNN 1 SNNCIOHNN 
ee 
NOPFAOANMTOORNMSTNO ow 
i i oo oo oe | 


FFPFTOMNWINAYO00D000 UF 9D 


ee ee ee Ne a ee ee ee es ee ee ee tee ee <LI FEL 


aadadddadaecdadstoew 


? 
( 


A( 1,6) 

C 

@ 

5 
Dg) 
M6, 


Sololelololololelelelelejleleleloleleloleleleleleole) 
ANOFINOM DHOANMSFOOE DHDOAAUMSTINO 
OWDODODODOD0O0REE EEE erroonnnwmea 
Set st Ht tte 
lolol colelelcolele)|eolelelejeleloleleleleloelele) 


N ~ 
=~ es = 
= Co Co 

en N ay 
oO} et ~ 
= © ~ < 
ears rt m3 
ol et N= O 
zl - mJ ee 
Sa ~ ° ie 
=_—— = N ie) +t4 <a 
Cos =) i t+ O~ Xx 
wk « — oO - . | 
vin N a) Oo Z2ZNNZ= 
oOat~— ~— « « st “N= a> 
wt 2 zie (\aje=s, (a e FA Oe oe 
Wo owe ese a Am~ re 1X1.Y 
=O i) wot Wirt i NW ad th ow 


RRO OFYMISe DADVIA a9AY eM D 
DUMO eay (Qa FeSs re IOZ+t Mil ~Za 
OCOe2Z et ONIMO HS ITYI7D eo] It rONte 
CIWANO etNIY~ eA «Q J IOI oF 
ona= it Oo ~OlZwe LoS ied Fr YS 
DZ TOO WOW WHAOWO Ilnmn OL COU AO~O 
N= QOOOYNORSNYORONTATORIORYOTO 


(e) ot IN CS) + Ce N MOS 
ret oo aA es a NM mo mM 


60 


’ 


fi - t 
tPYY f2o5cef 
(CUT 20IG ef ine 


lelelelele) 
~-Oonoe 
DWDDONM 
Sa et 
lejelelele) 
oCo00o 


61 


> 
’ 


APPENDIX E 
PROGRAM VIPIPE 


A. GENERAL 

As a part of the present thesis investigation, attention 
was given to a digital computer program called VIPIPE which 
was originally written by Fink [1] and later modified and 
augmented by Kim [2]. This program provides solutions for 
in-plane and out-of-plane vibration of a single plane piping 
configuration having quite general geometry (except, of 
eourse, for its limitation to a single plane). It was 
desired to convert this program from a FORTRAN 63 program 
for use with CDC-1604 computer formerly employed at the 
Naval Postgraduate School, to a FORTRAN IV program for the 
IBM 360/67 presently in use, in order to provide an indepen- 
dent means of checking the results obtained as indicated in 
the body of this thesis and, also, because it would be 
useful to have a program of such generality available. For 
G@ecails of the theory of VIPIPE, refer to the theses written 
Dy fank and Kim. 

Because of the limitations on word length in IBM 360/67, 
the converted program, which was called VIPIPE-l, has been 
written to employ double precision arithmetic. Certain 
improvements have been made in the capability to produce 


graphical output. 
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The effort to make this conversion has been found to be 
much greater than was anticipated. It seems that CDC-1604 
computer is a more "forgiving" machine than is the 
IBM 360/67. In checking out the new program certain bugs 
came to light which were actually present in the earlier 
Version. but which apparently, did not..cause trouble... However, 
they represented some syntactical errors that resulted in 
failure on the IBM 360/67. In some cases it was clear how 
Go modify the. program so. as to eliminate, the. errors bul. in 
@ther cases the difficulty was so deep in the logic that a 
correction could not be made without acquiring as deep a 
knowledge of the program as was possessed by original 
@racers, Fink and Kim. 

Thus it cannot be asserted that VIPIPE-1 has been 
Ssuecessfully converted in every fine detail and it is 
recommended that additional efforts be made to assure the 
overall integrity of the program and to assess its limita- 
Zoos. Lt should be noted, however, that no gross errors 
have been found in any results obtained by use of the new 
peoeram:; the only question is that of the degree of accuracy 
Bian At 1S possible to obtain. 

Coliparisons Of accuracy.are Gifficult to make.) In,the 
forests place ~ the word-length is different in the new program, 
and, other things being equal, one would expect the new 
program to provide better accuracy , and t.o, maintain. ats 
integrity in the face of round-off error beyond the point 


where the earlier programs "bombed out." Thus, agreement 
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between Kim and Fink and disagreement with VIPIPE-1 does not 
mean that the former are "correct" and the latter is in 
error. Furthermore, Kim reports having discovered and cor- 
rected some bugs which he found in Fink's version. Also, 
there is some difficulty in understanding the precise 
details of some of the problems treated by Fink and Kim. 
With this preface which admonishes the reader not to 
make superficial inferences from comparisons, in Section C 
some comparisons between results obtained by Kim and results 


obtained using VIPIPE are presented. 


B. DATA DECK ASSEMBLY 
The order of data deck is 
i-ieGraph titde data cards . 
ej orropdem statement) card: 
3. Boundary condition data card(s). 
4, Branch joining angle data card(s). 
5. Component extensive properties data card(s). 
6. Component intensive properties data card(s). 
Field specification and format of graph title cards are 
given in Section E-l-c of Appendix B and information for the 


remaining cards are given in [1]. 


C. SAMPLE PROBLEMS 
Properties of piping systems analyzed in this section 
are given in [1]. Here, only numerical values of frequencies 


will be given. 
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1. System 1 (in-plane) 


Mode Comparison®* WiPibeo WEPLPRoT 
No. Values Delta m. (D) Delta m. 

A 75-09994540 75.09994544 7509994535 
2 207.01589137 207.01589036 207.01589099 
3 405 .63391916. 405. 83393183 40583391910 
4 670. 86408383 670. 86310856 670. 86408365 
5 1002.15520225 1002.19978814 1002.15520252 
6 1399. 70436269 1399. 70434917 
7 1863.51172599 1863.51198134 


Sources of comparison frequencies for system 1 and system 2 
are piven in Section D. 


*#*#Frequencies are given in [1]. 


#** Frequencies are found by VIPIPE-1. 


2. System 2 (in-plane) 


Mode Comparison VIPLPE VIPIPE-1 
No. Values Delta m. (D) _ Delta m. 
a) 11. 80213586 11. 80213389 11. 80213432 
2 73.96272297 13296272007 73.96272087 
3 207.09776595 207.09776685 207.09776559 
4 405.82896584 405.82894132 405.82896568 
5 670.86435913 670.86200817 670.86435889 
6 1002.15518770 1002. 26833752 1002.15518758 
7 1399. 70436229 1399.08849534 1399. 70437163 
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3. System 4 (in-plane) 


Mode 
No. 


ine) 


—- W 


1S] 


VIPIPE 
M-D m. (D) 


66.699 34902 
171. 22352538 
335. 35007063 
535.64529414 
773. 44358487 


4, System 5 (in-plane) 


Mode 
No. 


1 
2 
S 


Mode 


NN nN WW LF WW 


VIPIPE 
Delta m. (D) 


116.89479180 
304.88354398 
566.43670976 


VIPIPE 
BaM # my 1D) 


826. 88281929 


1126.83797508 


1666.58878192 
2099.52844030 
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VIPIPE-1 
M-D m. 


67. 31642169 
172. 86857155 
339 .09491357 
544.08929568 
794 .57343086 


VIPIPE-1 
Delta m. 


118.95185774 
302.91288556 
519). 19656352 


VIPIPE-1 
Bais im, we 


118.95185779 
302.91288592 
519.19656389 
794.49003694 

1126 .27858174 

1464.59530724 

1874.97923215 


D. SOURCES OF COMPARISON FREQUENCIES FOR SYSTEM 1 AND 
SYSTEM 2 


The governing fourth order equation of a homogeneous 


Single component straight system 
q ; 
Ye (E-1) 
was solved to give 


y = aCosPx + b Cosh Px +-¢ Sin Px-t d Sinh Px 


Upon applying boundary conditions for system 1* Grey = 0 
ioexe= O.. X = 1) and simplifying, 16 wes found that the 
eigenvalues of the system must satisfy 


Cos(€PL.) Cosh(eL) = 1 =) 


Upon applying the boundary conditions for system 2 (y=y'= 0 
pee = 0, y" = yy" = 0 at x = L), it was, found that the 


eigenvalues of the system must satisfy 
Cos({PE)) Cosh(PL) = =1 (E-3) 


The roots of equation E-2 and E-3, from which the natural 
frequencies of corresponding systems may easily be found, 


were obtained through iteration using a digital computer. 


*System 1 is clamped-clamped end beam 
System 2 is clamped-free end beam 


67 


BH. PROGRAM LESTING 
il. inde 


Main program: 
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Section Sequence number 
INPUT 000010 -— 001080 
INVARIANTS 001120 - 001520 
CONTROL IN PLANE 001560 - 002470 
CONTROL OUT OF PLANE 002520 - 003420 
ITERATION 003460 - oo4460 
OUTPUT oo4490 - 008135 
Subroutines: 
SUBSEC 008140 -— 008640 
DISTM 008670 - 009220 
DISTMO 009652 - 009676 
RIGID 009652 - 009676 
RIGIO 009679 - 009730 
STIFCO 009760 - 009970 
STIFOO 010000 - 0102830 
STAVEC. 010260 -— 010550 
STAVEO 010580 - 010830 
INVERT 010860 =—_011570 
HANGER 011600 = 011770 
HANGEO 011800 - 011970 
CFIELO 012000 - 012520 
POINO O22550 = O127 76 
CFIELD 012000 - 012520 
POINT 012550 - 013490 


es 


SFIELD 


SFIELO 


MATMUL 


FINBRA 


FINMAT 


BRANCH 


BRANCO 


STATEM 


STATEB 


BRCOR 


DELMA 


Listing 
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023520 
013910 
014210 
014360 
014650 
014950 
015610 
016320 
016870 
017630 
017970 


013880 
014180 
014330 
014620 
014920 
015580 
016290 
016840 
017610 
017940 
019380 


100 
101 
Oak 
OZ 
102 
103 
104 
106 
Ont 
108 


109 
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